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Abstract — The performance of the V-BLAST approach, which 
utilizes successive interference cancellation (SIC) with optimal 
ordering, over independent Nakagami-wi fading channels is 
studied. Systems with two transmit and n receive antennas are 
employed whereas the potential erroneous decision of SIC is also 
considered. In particular, tight closed-form bound expressions are 
derived in terms of the average symbol error rate (ASER) and 
the outage probability, in case of binary and rectangular A^-ary 
constellation alphabets. The mathematical analysis is accompa- 
nied with selected performance evaluation and numerical results, 
which demonstrate the usefulness of the proposed approach. 

Index Terms — Nakagami Fading, V-BLAST, Successive In- 
terference Cancellation (SIC), Multiple Input-Multiple Output 
(MIMO), Bit Error Rate (BER), Correlated Fading. 



I. Introduction 

THE V-BLAST approach represents a cornerstone recep- 
tion strategy for multiple input-multiple output (MIMO) 
infrastructures because it achieves a high spectral efficiency 
and a substantial capacity gain IT], El- It utilizes successive 
interference cancellation (SIC) in a number of consecutive 
stages. The symbol detection and the corresponding decoding 
at a given SIC stage can be implemented according to an 
optimal symbol ordering, based on the highest signal-to-noise 
ratio (SNR) level, or without ordering. Since SIC is quite 
a complex process, its average symbol error rate (ASER) 
performance has been studied mainly numerically (e.g. Monte 
Carlo simulations) and/or semi-analytically with respect to the 
instantaneous symbol error rate (SER). 

Thereby, analytical research studies for the V-BLAST (or 
SIC) approach are very limited in the bibliography so far. More 
specifically, Loyka et al performed an analytical framework 
with respect to ASER for 2 x n MIMO systems with optimal 
ordering in [3| and for the generalized I x n case without 
optimal ordering in [2], where I and n denote the number 
of transmit and receive antennas, respectively. Nevertheless, 
these contributions assumed an error-free SIC approach and 
Rayleigh channel fading conditions. 

Nakagami-m is a versatile model, which includes the 
Rayleigh fading condition as a special case. To this end, an 
analytical framework for 2 x n MIMO SIC-enabled systems 
with optimal ordering over Nakagami-m fading channels is 
presented into this letter The merits of the proposed ap- 
proach are twofold: 1) tight closed-form bound formulae for 
V-BLAST systems over spatially independent Nakagami-m 
fading channels are derived in terms of ASER and the outage 
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probability, thereby generalize some of the results given in 
iH, 0; 2) a novel analytical expression for the potential error 
propagation of the SIC process is presented. 

II. Statistics of the SIC Stages 

Consider a 2 x n MIMO SIC-enabled system with two 
transmit and n > 2 receive antennas. The following standard 
baseband discrete-time system model is employed, expressed 
as r = H s + w, where H = [hi,h2] denotes the nx2 channel 
matrix, r = [ri, r™]"^, s = [si,S2]^, w = [wi, ...,Wn]'^ 
are the received, the transmit and the additive white Gaussian 
noise vector, respectively. Moreover, hj represents the ith 
channel nxl column vector, i = 1, 2 and (.)^ denotes vector 
transposition. 

Let X be the received instantaneous SNR. The probabil- 
ity density function (PDF) and the cumulative distribution 
function (CDF) of x over spatially independent Nakagami- 
m fading channels are, respectively, expressed as fx{x) = 
i^^M^x™-! sxp{-^x) and T,{x) = 7 (to, f a:) /r(m), 
where r(.) is the gamma function |5j eq. 8.310/1], 7(., .) 
is the lower incomplete gamma function ||5] eq. 8.350/1], 
TO = 2nxmN is the normalized Nakagami-m fading parameter 
with respect to the number of transmit and receive antennas 
(the factor 2 indicates the number of transmit antennas), 
rriN > ^ is the distribution shape parameter (which indicates 
the fading severity), Q, — £[x] is the average signal power and 
£[] denotes expectation. 

The reception process is implemented successively in a 
number of stages, proportional to the number of the transmit 
antennas. The key idea is that as a given symbol is detected, 
decoded and then canceled from the composite signal at the ith 
SIC stage, the remaining signal at the next stage experiences 
better channel conditions in terms of SNR and, hence, a better 
ASER performance. For a detailed description of the SIC 
architecture and methodology, see [[61. 

A. First SIC Stage 

The upper bound on the CDF of x given at the first SIC 
stage, Fi{x), is obtained as |i4( eq. 5] 

Fiix) - 



[n - 1) T.. (I) e-'dt >{n-l) Tl (|) f'-^dt, 



(1) 

where the expression on the right hand side of the inequality 
represents the actual CDF of x at the first SIC stage, which 
is analytically infeasible to be obtained in a straightforward 
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closed-form solution, mainly due to the involvement of 7(., .)^. 
However, unlike the actual CDF of x, the respective upper 
bound can be derived in a closed-form expression. Based on 
([B while utilizing first eq. 9.31/2] and then E eq. 7.81 1/2] 
we have that 
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where G™^" 



. ] is the Meijer's G function 15] eq. 9.30]. 
Thereby, taking the first derivative of (|2), the PDF of x 
given at the first SIC stage is expressed as 

dFi{x) 



/i(^) 



dx 

r(m) 



TO 
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where r(.,.) denotes the upper incomplete gamma function 
IS eq. 8.350/2]. 

B. Second SIC Stage 

The CDF of x given at the second SIC stage, F2{x), is 
obtained as 

F2(a;) =min[a;i,a;2] = 1 - [1 - T,{x)f = 2 T,{x) - T^ix). 

(4) 

In fact, the outage probability at the second stage is twice the 
corresponding one at the first stage [3, eq. 35], i.e. 2J^r(2a:^) — 
J-'^{2x). This occurs because the post-processing noise power 
at the second stage is twice of the branch noise power level, i.e 
af = {n — 2 + i) a^, where af and denote the noise power 
at the /th SIC stage and the total noise power, respectively. 
For a detailed description of this effect, see |[3] Appendix I]. 
Thus, we have that 

F2{x) « 2 T^{2x) = 7 f m, ) , ^ cx). (5) 

1 (to) \ It J 

It should be noted that the approximation of (|5]l converges to 
the actual ^2(0;), as given in (|4|, for typically medium/high 
average SNR values. However, the above mentioned fluctu- 
ation is maintained small even in the low SNR regime, as 
demonstrated by the numerical results of the next section. 
Hence, the PDF of x given at the second SIC stage is derived 
as 

om+l { m\^^ 
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III. Performance Analysis 

A. ASER 

The average symbol error rate (ASER) at the /th SIC stage, 
Ps.i, can directly be evaluated by averaging the conditional 
error probability (CEP), Pe{e\xi), over fi{xi) and is expressed 

as _ 00 

Ps.i = / Pei^lXt) fi{Xi)dXi. (7) 



(8) 



In case of binary modulations, CEP is defined as Q 

r(/3, axj) 



Peie\xi 



where a and /3 are certain constants that define the modulation 
type. In case of rectangular A^-ary modulations and higher 
values of the average input SNR, CEP is defined as [TJ 



P„(e\xj) = a erfc 



(9) 



where erfc(.) denotes the complementary error function [5, eq. 
8.250/4]. 

The total ASER in a 2 x n MIMO SIC-enabled system is 
expressed as 



Ps.total = Ps.i + P.s,2i^ - Ps.i) - Ps.i + Ps,2 " Pe,l * Pe,2. 

(10) 

It is worth noting that ( fTOl i is rigorous and accounts for the 
potential error propagation of the SIC process (i.e when no 
error occurs at the first stage and an error occurs at the second 
stage or both the SIC stages are erroneous). Unlike the first 
two terms of ( fTOl l, the third term involves a conditioning on 
both xi and X2, which are not statistically independent. Hence, 
the more complicated bivariate (correlated) PDF is required 
in this case. Moreover, P^.i * Pe,2 represents a second order 
statistic, which may fluctuate the performance of the total 
ASER, especially in the low SNR regions (where the presence 
of the error propagation is more emphatic). In the following, 
tight closed formulae are derived with respect to the total 
ASER, for binary and A^-ary modulation scheme^. 

1 ) ASER at the First SIC Stage: Based on Q and ^ while 
utilizing IB] eq. 2.10.6/1], the ASER for binary modulations is 
derived, as given in ( fTTT i. where pFq denotes the generalized 
hypergeometric function |5, eq. 9.14/1]. Similarly, based on 
Q and (|9]l, the ASER for A^-ary modulations is obtained in 
([T2]l, by invoking |8i eq. 2.10.8/2]. 

2) ASER at the Second SIC Stage: In case of binary mod- 
ulations, based on (|6]l and (|8]l while utilizing [8] eq. 2.10.3/2], 
the corresponding ASER is obtained, as given in (fT3T l. In case 
of A^-ary modulations, based on (|6]l and (|9]l while invoking 
|8, eq. 2.8.5/6] and after performing some straightforward 
algebraic manipulations, the corresponding ASER is derived 
in (O. 

3) Cross-Product (Correlated) ASER: The statistically cor- 
related cross-product term can be obtained by averaging CEP 
over the PDF of such an event and is defined as 



,l*Pe,2^ I 

Jo 



^e(e|?/) fyiy)dy, 



(15) 



where y = X1X2, Pe{e\y) denotes the CEP on y and fy{y) is 
the cross-product PDF, since (|5]l and the following condition 
hold [9 eq. 6.74] 



fyiy) 



/ fxi 
Jo XI 



(16) 



2r(/3) ' 



'Note that the average bit en'or rate at the (th stage, P5 i, can be easily 
deduced from the corresponding ASER, assuming that Pi, i = P^ i/logj A1. 
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The bivariate Nakagami-m 
as eq. 2] 



j(a;i,a;2), is expressed 



4(a;iX2) 



X I, 



r(m)f]™+i(l - p)p- 



(17) 



where denotes the modified Bessel function of the 

first kind and order v f8^ Appendix II. 10/1] and p repre- 
sents the correlation coefficient, which is denoted as p = 
cov(a;i, xi)l •\/var(xi)var(a;2), where var(.) and cov(.,.) de- 
note variance and covariance, respectively. By invoking [5 
eq. 3.478/4], ( fT6b can be easily resolved as 



Jy{v) 



X X 



4'"r(m)f^'"+i(l - p)p^ 

py 



JC 



(18) 



where IC^{.) denotes the modified Bessel function of the 
second kind and order i' [8 eq. II. 10/2]. 

Unfortunately, substituting (fTsT l and (O or (|9) in ([15) does 
not provide a straightforward tabulated closed-form solution. 
Nonetheless, by applying the infinite series representation of 
||8j Appendix II. 10/1] in ^TE^, while utilizing |8, eq. 
2.16.1/31], the correlated ASER for binary modulations is 
obtained in ( fT9] l. where denotes the digamma function |5 

^It is assumed that Qi = Qi = f^. the pre-processing average SNR is 
identical for all branches and for both transmitters. 



eq. 8.360/1]. In addition, by using ^ and ( fTSl l in ( fTSb and with 
the aid of O eq. 2.16.1/25], the corresponding ASER for M- 
ary modulations is derived in ( |20l l. Note that the summation 
terms in (fT9] l and ( |20] i converge very rapidly for various system 
scenarios (e.g. (l20l l requires only 4 summation terms in order 
to converge up to the 9* digit, when m = 2 and p — 0.7). 

Figs. [T] and |2] show the ASER performance of the two 
consecutive SIC stages and the total ASER, respectively, in 
various system scenarios. It is obvious that the corresponding 
ASER is sharply affected with an increase of the number 
of receive antennas and with a reduction of the channel 
fading severity (e.g. in higher mjv values). Moreover, in the 
worst system scenario in terms of ASER, where only two 
receive antennas are employed and the channel fading severity 
is quite intense (mpf = 0.5), the influence of the cross- 
product ASER is shown in Fig. [3] The numerical results at the 
above mentioned configurations have obtained via numerical 
evaluation based on the actual statistics, as given at the right 
hand side of ([TJ and dUi, and then by utilizing ( fTob . A slight 
difference on the ASER performance between the analytical 
bound formulations and the respective exact numerical veri- 
fication is observed whereas quite an effective computational 
gain is achieved by performing the proposed approach. 

B. Outage Probability 

The outage probability, Pout,i{xth), at the first SIC stage is 
directly obtained from ([2), where xth denotes a threshold SNR 
value. The outage probability at the second stage, Pout.zixth), 
conditioned on an error-free first stage is given in (2}. The 
corresponding unconditional outage probability at the second 
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SNR (dB) 

Fig. 1. The ASER of the first and the second SIC stage vs. different values 
of the average input SNR, for a 2 X n MIMO system, when mjv = 1 and 
a BPSK modulation scheme is considered. The numerical and the analytical 
results are indicated by marker signs and lines, respectively. 
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Fig. 4. The outage probability vs. different values of the average input SNR 
for a 2 X n MIMO system, when mjv = 2 and a BPSK modulation scheme is 
considered. The numerical and the analytical results are indicated by marker 
signs and lines, respectively. 
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Fig. 2. The total ASER, as given in OOK vs. different modulation schemes 
and several values of the average input SNR, for a 2 X n MIMO system, when 
mjv = 2. The numeiical and the analytical results are indicated by marker 
signs and solid lines, respectively. 
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Fig. 3. The cross-product ASER vs. different values of the average input 
SNR, for a 2 X 2 MIMO system, under various correlation conditions and 
different modulation schemes, when mjy = 0.5. 

Stage, Pout li^th), which considers the potential erroneous 
decision at the first stage, is expressed as Pgutzi^th) = 
F2{xth){i — Ps,i)+Ps.i- Note, that Pout,i is independent of the 
error propagation whereas Pl^^f_ 2 is typically upper bounded by 
Fig. ^indicates the outage probability in various average 
SNR regions for 2 x n MIMO systems. It is obvious that as 
the spatial diversity gain increases, the outage performance 
improves (i.e. Pout decreases) for both SIC stages. 

IV. Performance Bounds for the Generalized 
V-BLAST 

Consider an I x n MIMO system with I < n transmit 
antennas. Hence, H = [hi,...,h;], r = [ri, r„]^, s ~ 
[si, ...,si]'^, w = [wi, ...,Wn]'^ while SIC is implemented in 
I consecutive stages. In this case, the distribution of SNR can 
not be resolved in a closed-form expression for the ith stage. 



when i > 1 (even if the signal undergoes Rayleigh fading fTT] 
eq. 25]). Fortunately, there is an upper bound expression for 
the outage of the first SIC stage, given as ID 

fI'^-\x) = - 1) [.F. (f )]' - ty-^dt, 

Jo 

(21) 

which represents a generalization of ([T). Note that (I2II can 
not be evaluated in a closed formulation mainly due to the 
involvement of 7(., .) within Fx{x). However, comparing ([T]i 
and ( I2TI) whereas recognizing that the diversity gain of an Z x n 
V-BLAST at the first stage (which is n — Z + 1) is always lower 
than a 2 X n V-BLAST, i.e. n-l + 1 < n-1 for I > 2, we have 
that f/'^"' > Fi(('-i)><") > ... > f/'^"^ while equiprobably 
it holds that Pi|i^") > P«j'-i)><") > ... > pjf Thereby, 
(|2]l and (fTTl l (or ilT^ ) can serve as sharp closed-form lower 
bounds for the generalized I x n case with respect to the outage 
performance and ASER, respectively. 
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